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MAGNETIC  RECONNECTION  IN  A  NON-MAXWELLIAN  NEUTRAL  SHEET 


I.  Introduction 


The  collisionless  tearing  mode  (Furth,  1962;  Pfirsch,  1962;  Laval  et 
al,  1966)  has  attracted  considerable  attention  in  connection  with  the 
earth's  magnetotail  (Coppi  et  al,  1966;  Schindler,  1966).  It  was  suggested 
(Coppi  et  al ,  1966)  as  a  possible  mechanism  for  magnetic  field  reconnection 
in  the  magnetosphere.  The  tearing  mode  may  also  be  relevant  to  the  day3ide 
magnetopause  (e.g.,  Greenly  and  Scnnerup,  1981;  Guest  and  Coroniti, 
1981).  The  basic  idea  is  that  the  magnetic  field  energy  stored  in  the 
magnetotail  as  a  result  of  interaction  with  the  solar  wind  is  released  in 
other  forms  of  energy  via  reconnection  processes.  The  collisionless 
tearing  mode  is  a  possible  instability  that  allows  necessary  changes  in  the 
magnetic  topology  in  the  absence  of  resistivity. 


The  geometry  most  often  used  to  model  the  magr.etotail  is  the  neutral 

sheet  geometry.  The  standard  analysis  (e.g.,  Laval  et  al,  1966)  uses  a 

Maxwellian  plasma  distribution  function.  However,  in  a  collisionles3 

plasma  such  as  the  magr.etotail,  the  motion  of  particles  parallel  to  the 

magnetic  field  is  decoupled  from  the  perpendicular  motion  and  r.on- 

Maxwelliar.  features  are  likely  to  persist.  Laval  and  Pellat  (1  963)  shewed 

that  the  collisionless  tearing  mode  ( k  j  )  can  be  strongly  modified  by  weak 

electron  temperature  ar.isotrooy,  il  /T  -  1  i  <  p  /6  <<  1  ,  where  k  is  the 

e±  8 1  1  e  — 

wave  vector,  J3q,  i3  the  equilibrium  magnetic  field,  og  is  the  typical 

electron  Larmor  radius,  5  is  the  character i3tic  half-width  of  the  neutral 

sheet,  and  T  and  T  are  the  electron  temoeratures  associated  with  the 
ej_  ej 

particle  motion  perpendicular  and  parallel  to  3-,  respectively. 
Sucsequer.tly ,  Forslur.d  (’  963  )  obtained  an  approximate  dispersion  -elation 
with  weakly  anisotropic  electrons.  The  work  showed  a  substantial 

consistent  with 


enhancement  of  the  tearing-mode  growth  rate  for  T  > 

U 

the  result  of  Laval  and  Pellat  (1963). 
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In  the  werk  of  Forslur.d  (*968), 
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axis-crossing  icr. 
:  of  the  order 
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‘giigitie  (the  convent  i  oral  ''two-region" 
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approximation) .  Recently,  Chen  and  Paimadesso  (1984)  have  shown  the 
existence  of  an  ion  intermediate  region  with  a  thickness  of  the  order  of 
(p.5)1/2  >>  (pgo)1/2  where  the  axis-crossing  ion  orbits  make  a  major 
contribution.  Inclusion  of  the  ion  intermediate  region  (the  "three-region" 
approximation)  shows  that  a  given  degree  of  ion  anisotropy  (7^/7^  >  1  ) 
car.  increase  the  growth  rate  by  nearly  one  order  of  magnitude  over  the 
result  obtained  using  the  two  region  approximation.  Furthermore,  in  the 
twc-regicn  approximation,  because  the  ion  anisotropy  effects  are  neglected, 
a  neutral  3heet  with  f  1  —  T  /7  ]•  p  /5  «  1  is  completely  stable  unless 
unrealistically  high  degrees  of  ion  anisotropy  is  included.  However,  the 
three-region  treatment  shows  that  only  modest  ion  anisotropy  is  needed  to 
destabilize  such  a  neutral  sheet.  This  is  of  particular  interest  since  any 
significant  anisotropy  in  the  electron  population  would  tend  to  be 
isotropized  on  very  fast  time  scales  (see,  for  example,  Coppi  and 
Resenbluth,  1968).  The  time  scale  for  the  ions  is  longer,  of  the  order  of 
a  minute,  making  it  more  relevant  for  magnetotail  reconnection  processes. 

The  three-region  approximation  also  shews  that  a  more  adequate 
treatment  of  physically  interesting  systems  containing  anisotropic  ions 
requires  an  accurate  treatment  of  large  ion  orbits.  In  the  above  works, 
the  complicated  equilibrium  orbits  and  orbit  integrals  were  approximated. 
The  inability  to  treat  orbits  exactly  leads  to  a  number  of  limitations  such 
as  the  oonstar.t-’y  approximation  (i.e.,  1  is  constant  in  the  inner  region) 
and  weak  anisotropy.  In  addition,  if  anisotropic  ion  effects  dominate  in 
t .he  ion  intermediate  region  wnile  the  electrons  still  dominate  in  the  inner 
region,  then  the  orbits  for  both  species  with  disparate  scale  sizes  must  te 
creased  accurately.  Hoidrer.  (1970)  carried  out  full  integro-dif ferer.tial 
calculations  of  the  tearing  mode  in  a  (relativistic )  neutral  sheet.  In 
tni3  work,  the  oroits  were  calculated  numerically  ar.c  an  Iterative  method 
was  iseo  with  discrete  variables.  As  a  result,  the  numerics  required  were 


substantia- . 


'gr  qr g 


bevel oced  an 


method  which  treats  all  the  eauilitrium  oroits  exactly 


analytically  for  both  species.  This  method  maxes  it  possible  to  remove  tr.e 


above  ..mi 


ir.ey  prevized  an  accurate  treatment  of  nig. hi y  oc.r- 


?3.xwell  i  an  neutral  sr.eet  plasmas  with  minimal 


Pi 


VvV'.VV 


In  this  paper,  we  incorporate  the  recent  results  into  a  model  of 
magnetic  reconnection  in  the  magnetotail.  We  will  first  review  the 
theoretical  results  described  above  pertaining  to  the  non-Maxwell i an 
collisionless  tearing-mode  instability  in  a  neutral  sheet. 

II.  The  Non-Maxwell i an  Tearing  Mode  (NMTM) 


Consider  a  neutral  sheet  whose  magnetic  field  profile  is  given 


by  3^ ( z )  «  3x(z)  *  3uch  that  B^f-z)  “  ~  3x(z).  current  density  is 


given  by  J(z)  *  Jq(z)_£  »  (c/Uir)7  x  3^  with  (-2)  -  Jq(z).  The  equations 

of  motion  admit  three  independent  constants  of  motion;  H  »  (i/2)mv 

0  ±2  2  -L  2 

P  »  mv  +  (q/c)A  (z),  and  H  «  (1/2)  m  v  where  v  -  v  +  v 

yy  y  -l  I  ±  y  2 

and  v  »  v  .  Here  A(z)  »  A  (z)  y  is  the  equilibrium  vector  potential  and 
I  x  -  y  i 

the  equilibrium  electric  field  is  set  equal  to  zero.  The  tearing 


instability  oar.  be  described  by  perturbations  of  the  form 
•J»(x,s,t)  3  4»(z)exp(  ikx  -  iat)  where  the  wavevector  k  ■  kx  is  parallel  to 

the  equilibrium  magnetic  field.  Although  not  necessary,  we  will  consider 
perturbations  with  |u|  <  a  .  and  neglect  perturbed  scalar  potential.  We 
assume  charge  neutrality  to  first  order. 


We  consider 
of  the  type  ?,  * 
standard  method 
function 


a  class  of  equilibria  described 
?.(H  ,  ?  ,  H  )  where  j  is  the 

j  1  y  1 

of  characteristics,  the  first 


by  distribution  functions 
species  index.  Using  the 
order  Vlasov  distribution 


for  each  species  can  be  written  as 


The  integration  is  carried  out  along  equilibrium  orbits 
vector  potential  then  satisfies  the  equation 


The  perturbed 


ill 

c 
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<-•  / 
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C) 


where 


iy 


z 

j 
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J* 


As  a  general  remark,  the  full  solution  of  the  first-order  problem  oar. 
be  obtained  by  solving  the  integro-dif f erential  equation  (3).  In  the  work 
of  Forslund  (1968),  an  approximate  dispersion  relation  for  weakly 
anisotropic  electrons  was  obtained  using  (two)  velocity  moments  to  express 
the  orbit  integral.  Essentially  the  same  dispersion  relation  was  obtained 
by  Cher,  and  Palmadesso  (*985)  by  matching  the  inner  \z\  <  {2pa5)'/~  and 
outer  (|3|  >  (2po6)  *)  solutions  which  are  both  analytically  accessible. 
In  this  work,  the  straight-line  orbit  approximation  (Coppi  et  al . ,  1966; 
Cobrowoir.y,  ’968)  was  used  for  the  axis-crossing  orbits  in  the  electron 
inner  region  for  both  species.  Both  approaches,  however,  neglect  the  fact 
tnat  the  axis-crossing  ion  orbits  extend  far  beyond  the  electron  region. 
By  including  the  axis-crossing  ion  orbits,  Cher,  and  Palmadesso  '  *  985 ) 
snowec  tr.e  existence  of  an  ion  intermediate 


1  .'2 


,1/2 

(2p.6  , 

i 


•egion,  2  p  5 )  "  ■'  i  o  i  <  ,'23,6}  ",  in  which  the  anisotropic  ’.or. 
contribution  dominates.  By  matching  tr.e  solutions  in  tne  tnree  regions, 
they  four.c  that,  for  a  giver,  degree  of  anisotropy,  the  growth  rate  based  on 
the  three-region  approximation  oar.  exceed,  by  nearly  one  order  of 
magnitude,  the  growth  rates  obtained  by  the  two-region  approximation. 
Figure  *  summarizes  these  results,  showing  the  growth  rates  for  several 
values  or  .  with  isotrcoic  electrons.  ir.  particular,  curve 


f  T  'T.  with  isotrcoic  electror.3.  I. 

U  M 

*  .  5  which  is  obtained  usi-.g  tne  tnree-regi  on  approximation 


-1  1  I 

is  significantly  nig.-.er  than  Cu: 


1 


e  wr. icr.  .s  obtained 


.he  two- -eg. 


approximation  :  or 
Y  w  aval  er.  g  t  ns  an  a 
base . 


:e  same  va^ue  o:  .. 
'.siderablv  shorter  for 


Her 


t.-.at  the  maximum 


tnar.  tne 
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In  the  above  works,  various  approximations  were  used  to  treat  the 
orbit  integrals  such  as  the  straight-line  orbit  approximation,  constant-^ 
approximation  and  expansion  in  terms  of  velocity  moments.  These 
simplifications  constitute  severe  limitations  on  the  applicability  of  the 
analyses  to  physically  interesting  systems.  Chen  and  Lee  (1985)  developed 
an  integro-dif  f  erer.tial  equation  method  which  uses  all  the  equilibrium 
orbits  for  both  species  exactly  and  analytically.  They  showed  that  the 
dispersion  relation  for  |<u|  <  w  nas  the  general  form 


(4) 


wnere 


y  -  y/'ji  .  with  u  .  -  e3„/m.c  3_  and  k  -  k<5 .  Here,  D  is  a  universal 
ci  ci  0  1  0 


function  of  T  /?  ,  T.  /T.  and  T  /T.  .  Using  the  Galerkir.  method,  they 

-1  - 1  1  I  *1 

solved  the  integro-dif ferential  equation  (3)  and  obtained  the  the  universal 


function  D  and  the  eigenmode  structure.  Figure  2,  reproduced  from  Chen  and 


Lee  (’985',  shews  the  dispersion  relation  for  several  values  of  T  /7a^. 
The  preferred  wavelength  is  seen  to  be  reduced  significantly  although  the 


neglect  of  the  scalar  potential  did  not  allow  the  accurate  determination  of 
the  value  of  <  for  the  maximum  growtn  rates.  In  this  worx ,  both  species 
are  allowed  to  be  highly  nor.-Maxwellian.  As  a  result,  the  instability  is 
dominated  by  the  electrons.  A  nor.-Maxwellian  icr.  population  with  isotropic 
electrons  has  similar  behavior  but  with  smaller  enhancement  of  the  growth 
rate.  A  recent  numerical  simulation  stucy  (Ambroslar.o  et  al ,  ’986)  closely 
suocort3  the  behavior  described  above. 


III.  Magnetic  Reconnection  in  a  Nor.-Maxwellian  Plasma 


Equation  shews  that  the  tvoe  of  ncn-Maxwei 1  i an  features  must  be 


'l 


i 


v 


representable  by  distribution  functions  of  the  form  F  = 
since  the  first  two  terms  are  the  usual  tearing  terms.  Here,  C  is  a 

r.  the  above  discussion,  we 


constant  of  metier,  indecentiert  of  H  and  ?.,. 

1  7 

have  used  tr.e  case  with  C  =  H  .  Note  that  it  is  ocssible  to  construct  a 


distribution  of  the  f 


non-Maxwell i an 


orm 


However , 


tvoe 


r.on-Maxwelliar.  tearing  mode,  in  effect,  consists  of  two  parts;  the  usual 
tearing  mode  and  the  part  driver,  by  the  non-Maxwell  i  an  distribution  of 
particle  energy.  The  latter  process  is  generally  much  stronger.  This 
suggests  that  the  NMTM  can  rearrange  the  internal  particle  energy  as  well 
as  the  magnetic  energy  on  a  fast  time  scale  (possibly  as  high 
as  '(  -  !i)^.  )  .  The  instability  should  lead  to  the  formation  of  small  scale 


islands.  The  islands  then  should  coalesce  rapidly  (Finn  and  Kaw,  *977), 
which  is  a  well  documented  process  (see,  for  example,  Pritchet  and  Wu, 
'979;  3rur.el  et  al . ,  '982;  Tajiaa  and  Sakai ,  1  986 ) .  This  scenario  is  also 
supported  by  the  simulation  results  of  Aabrosiaro  et  al.  ('986).  In 


addition,  these  results  show  that  the  conversion  of  magnetic  energy  to 
kinetic  energy  takes  place  faster  and  more  efficiently.  We  note  that  the 
above  process  may  play  a  role  in  "forced  reconnection"  since  the  driving 
forces  may  preferentially  increase  the  perpendicular  temperature.  The 
magnitude  of  anisotropy  and  hence  the  growth  rate  depends  on  the  strength 
of  t.ne  forcing  function. 


We  now  consider  a  possible  role  the  non-Maxwell i an  tearing  mode  may 
play  in  a  physical  system  such  as  the  magnetotail.  Consider  a  Maxwellian 
equilibrium  neutral  sheet.  Suppose  the  particle  distribution  function 

undergoes  some  large-scale  changes  due  to  cnar.ges  in  the  external 
cor.ditior.3.  This  leads  to  a  rearrangement  of  particle  energy.  If  the 
system  admits  3  independent  invariants  of  motion  as  cescriced  in  the 
precsecmg  paragrapns,  then  the  free-er.ergy  associated  with,  the  non- 
maxwell. an  features  can  drive  the  WMTM .  This  occurs  along  wit n  t.ne  usual 

tearing  mode.  For  typical  magnetotail  parameters,  the  linear  e-foldmg 

time  car.  ce  -educed  tc  a  small  f -action  of  a  minute  depending  on  the  r.c.n- 

Maxwel.  .  ar.  features  Cher  and  Pal. mad  esse,  '  ?cJ ;  Cher,  and  lee,  '985'. 


Subsequently ,  3mall  scale  C<6  >->  *  )  islands  are  formed  wr.ich  tr.er.  coalesce 


.  e  •  ima  a 


:ly.  The  er.t.re  process  may  manifest  itself  as  a  triggering  of 
irection  following  c.har.ges  m  pr.ysical  conditions.  It  may  also  appear 
e  explosive  .  laleev  ei  al .  .  '}"2;  Terasawa,  '?3‘;  loror.ioi.  ’985: 

Saw. a ;  .  '  0P6  ;  . 


In  the  context  of  a  3ubstorm,  the  substorm  may  appear  to  be  triggered  or 
directly  driven  in  response  to  some  changes  in  the  solar  wind  conditions 
with  the  NMTM  providing  a  mechanism  for  the  "solar  wind  control".  The  time 
delay  is  determined  by  the  fast  non-Maxwell i an  tearing  time  scales  and  the 
subsequent  coalescence  of  islands.  It  is  interesting  to  note  that  there 
has  been  considerable  effort  to  classify  substorms  in  two  basic  categories; 
storage-unloading  of  energy  and  direct-drive  by  the  solar  wind  (Akasofu, 
1981)  The  above  mechanism  contains  attributes  of  both  types.  In  this 
regard,  the  differential  memory  of  particle  distributions  suggested  by  Chen 
and  Palmadessc  (1986)  provides  a  natural  means  for  the  generation  of 
necessary  ncn-Maxwellian  distributions. 

Finally,  it  has  been  suggested  (Galeev  and  Zeleyni,  1976;  Lembege  and 
Peliat,  1982)  the  magnetic  field  component  normal  to  the  equatorial  plane 
can  modify  the  tearing  mode  properties  significantly  if  the  electrons  are 
magnetized,  possibly  stabilizing  the  ion  tearing  mode  (Schindler,  1979). 
However,  Coror.iti  (1930)  has  pointed  out  that  pitch  angle  scattering  of 
electrons  car.  effectively  demagnetize  the  electrons  so  that  the  tearing 
mode  can  still  be  unstable.  Recently,  Chen  and  Palmadesso  (1985)  have 
shown  tr.at  the  normal  component  renders  the  single-particle  motion  in  a 
quasi-neutral  sheet  ncnintegrable  sc  that  the  electrons  may  be  only 
partially  magnetized.  In  addition,  Goldstein  and  Schindler  (1973)  and 
Swift  (1?35)  have  suggested  that  the  ionospheric  influences  may  also 
contribute  to  the  oestabilizati  cr.  of  the  ( conventional )  tearing  mode. 
Thus,  tr.ere  are  3tror.g  indications  that  the  nor.-Maxweilian  tearing  mode  is 
important  fcr  magnetotail  dynamics. 
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